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c/3 . Abstract 
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Ether-drift experiments are the only known experiments which, in principle, can distin- 



guish Einstein's interpretation of the relativistic effects from the Lorentzian point of view with 

some underlying form of ether which plays the role of a preferred reference frame. To test 

experimentally this latter point of view, it is usually assumed that the macroscopic Earth's 
> . 

, motion should be detectable in the laboratory from the time dependence of the data. There- 



fore any observed stochastic signal, which does not exhibit the smooth modulations expected 



■ from the Earth's rotation, tends to be considered as a spurious instrumental effect, e.g. ther- 

' mal noise. The real situation, however, might be more subtle if the hypothetical ether (i.e. 

^SJ ■ the physical vacuum) resembles a turbulent fluid where large-scale and small-scale motions 

are only indirectly related. In this case in fact, besides thermal noise, the data might contain a 
genuine physical stochastic component. To test this scenario, we have performed a numerical 
simulation to estimate the signal in a toy model of a statistically isotropic and homogeneous 
turbulent flow. After subtracting the known forms of disturbances, the present data are con- 
sistent with velocity fluctuations whose absolute scale is determined by the Earth's cosmic 
motion with respect to the CMB (projected in the plane of the interferometer at the lati- 
tude of the laboratory). Therefore the Earth's motion, although undetectable from the naive 
time-dependence of the data, could nevertheless show up in their statistical distributions. In 
particular, the predicted non-gaussian nature of the distribution of the instantaneous data 
will be precisely tested with the forthcoming generation of precise cryogenic experiments, 
with potentially important implications for our understanding of both gravity and relativity. 



1. Introduction 



Ether-drift experiments, starting from the original Michelson-Morley experiment, have played 
a fundamental role in the history of physics. Today, they are still crucial in view of the 
basic ambiguity which persists in the interpretation of the relativistic effects. In fact, the 
basic Lorentz transformations, rather than originating from the relative motion of a pair of 
observers S' and S", as in Einstein's relativity, might be associated with their individual 
velocity parameters /?' = v' /c and /?" = v" /c [H El [3] and thus reflect the existence of 
some underlying ethereal medium. Still, due to the fundamental group properties, the two 
frames S' and S" would also be mutually connected by a Lorentz transformation with relative 
velocity parameter 

^-1 - 1 _ i3>/3" - c 

(we restrict for simplicity to one-dimensional motions) . In spite of the deep conceptual differ- 
ence, this would produce a substantial quantitative equivalence with Einstein's formulation 
for most standard experimental tests, where one just compares the relative measurements of 
a pair of observers. Thus, to resolve the ambiguity, one should also consider the ether-drift 
experiments. 

At the same time, if the speed of light in the vacuum measured on the Earth's surface, 
say c^, coincides with the basic parameter c entering Lorentz transformations, relativistic 
effects conspire to make undetectable the individual velocity parameters /?', /?",... Therefore 
the only possibility is that c-^ and c do not coincide exactly^ see e.g. [Ij. In this case, in fact, 
there would be a tiny ether-drift effect 5 ~ /3^(c — c^)/c. We emphasize that this possibility, 
explored in ref . [5] in the framework of the so called emergent-gravity scenario [6l [7j , where 
(c — c^)/c was estimated [5] to be 0(10"^), is not ruled out by the present experimental data. 
In fact the ether-drift, as measured from the fractional beat signal between two vacuum 
optical resonators [8j, gives 8 ~ 10^"*^^ and thus could indicate an absolute Earth's velocity 
of about 300 km/s, as in most cosmic motions. This basic point can be easily checked by 
looking at Fig. 9(a) of ref. [9] where a typical sequence of 40 data collected at regular steps of 
1 second is reported 0. As one can see, this instantaneous signal exhibits random fluctuations 
of about ±1 Hz and this value, for a laser frequency of 2.82 • 10^^ Hz, might correspond to a 
genuine ether-drift effect of about ±3.5 • 10~^^. 

Still, so far, all experimental groups tend to interpret this type of stochastic signal as a 
mere instrumental effect, e.g. thermal noise. However, as discussed in [5], this interpretation 

^With respect to other articles, ref. [9] has the advantage to report the instantaneous raw data. The 
experiment also adopts a sophisticated geometrical set-up where, to minimize all possible asymmetries, the 
two optical cavities are obtained from the same block of ULE (Ultra Low Expansion) material. As such, the 
results of ref. ^ will play an important role in our analysis. 
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is not so obvious. For instance the same typical magnitude 10~^^ was obtained in [5] from 
an independent re-analysis of the cryogenic experiment of ref.|10). Since thermal noise is 
estimated on the basis of the fluctuation-dissipation theorem, there is no real reason that 
both room temperature and cryogenic experiments exhibit the same instrumental effects. 

To explore this important aspect, in this paper, we shall adopt the same point of view of 
ref.[5]. Namely, besides thermal noise (which is expected any way in all room-temperature 
experiments) there might be a genuine physical stochastic component in the data as, for 
instance, if the physical vacuum were similar to a form of highly turbulent ether. This idea 
is deep rooted in the basic foundational aspects of both quantum theory and relativity (see 
and finds additional motivations in those representations of the vacuum as a form of 
'space-time foam' which indeed resembles a turbulent fluid [l2l [TU El [JSj [l6]. For this 
reason, within a simple model of statistically isotropic and homogeneous turbulence, we have 
performed a numerical simulation to estimate the effects that could be detected in the next 
generation [IT] of precise cryogenic experiments. 

After this general introduction, the plan of the paper is as follows. In Sect. 2 we shall 
remind the problem of measuring the two-way velocity of light and in Sect. 3 the basics of the 
present ether-drift experiments. In Sect. 4 we shall describe our numerical simulation and give 
definite predictions for the next series of experiments. Finally Sect. 5 will contain a summary 
and our conclusions. 



2. The two-way velocity of light 

The basic concept in ether-drift experiments is the two-way velocity of light in the vacuum 
Cy{9). This is defined in terms of the one-way velocity c-y{9) (which is not unambiguously 
measurable) through the relation 

_ 2c^{9)c^{7T + 9) 



and could exhibit a non-zero anisotropy 

Ace _ c^{tt/2 + 9) -c^{9) 
c (c^) 



(3) 



This theoretical concept is related to the measurable frequency shift, i.e. the beat signal, Az^ 
of two optical resonators [8] through the relation 

where f o is the reference frequency of the two optical resonators and the suffix "phys" indicates 
a hypothetical physical part of the frequency shift after subtraction of all spurious effects. 
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As a possible theoretical framework for a non-zero anisotropy, we shall concentrate on a 
scenario which introduces some difference with respect to standard General Relativity and 
has a very simple motivation: c-y might differ from the basic parameter c entering Lorentz 
transformations due to gravitational effects. To this end, as anticipated in the Introduc- 
tion, one can consider the emergent-gravity scenario [6l [7] where the space-time curvature 
observed in a gravitational field becomes an effective phenomenon in flat space, analogously 
to a hydrodynamic description of moving fluids. In this perspective, gravity produces local 
modifications of the basic space-time units which are known, see e.g. |18tll9j. to represent an 
alternative way to introduce the concept of curvature [§. This scenario represents the simplest 
modification of the standard picture which allows for a non-vanishing anisotropy and gives 
the correct order of magnitude 6 ~ 10~^^. As such, it will be adopted in the rest of this 
paper. 

To better appreciate the possible implications for the measurement of the speed of light, 
it is convenient to start, as in ref. [5], from the basic notion: the definition of speed as 
(distance moved) /(time taken). To this end, one has to choose some standards of distance 
and time and different choices can give different answers. Therefore, we shall adopt the 
same point of view of special relativity: the right space-time units are those for which the 
two-way velocity of light in the vacuum c^, when measured in an inertial frame, coincides 
with the basic parameter c entering Lorentz transformations. However, inertial frames are 
just an idealization. Therefore the appropriate realization is to assume local standards of 
distance and time such that the identification c-y = c holds as an asymptotic relation in 
the physical conditions which are as close as possible to an inertial frame, i.e. in a freely 
falling frame (at least by restricting to a space-time region small enough that tidal effects 
of the external gravitational potential Ucxti-c) can be ignored). This is essential to obtain 
an operative definition of the otherwise unknown parameter c. With these premises, light 
propagation for an observer S' sitting on the Earth's surface can be described with increasing 
degrees of approximations [5l [H] : 

i) In a zeroth-order approximation, S' is considered a freely falling frame. This amounts 
to assume = c so that, given two events which, in terms of the local space-time units of 
S', differ by (dx , dy , dz , dt) , light propagation is described by the condition (ff= 'free-fall') 



^This point of view has been vividly represented by Thorne in one of liis books [20j : "Is space-time really 
curved ? Isn't conceivable that space-time is actually flat, but clocks and rulers with which we measure it, 
and which we regard as perfect, are actually rubbery ? Might not even the most perfect of clocks slow down 
or speed up and the most perfect of rulers shrink or expand, as we move them from point to point and change 
their orientations ? Would not such distortions of our clocks and rulers make a truly flat space-time appear 
to be curved ? Yes". 




(5) 
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ii) However, is really the Earth a freely-falling frame ? To a closer look, in fact, an 
observer S' placed on the Earth's surface can only be considered a freely-falling frame up to 
the presence of the Earth's gravitational field. Its inclusion leads to tiny deviations from the 
standard Eq.dS}. These can be estimated by considering S' as a freely- falling frame (in the 
same external gravitational field described by Ue:,^tix)) that however is also carrying on board 
a heavy object of mass M (the Earth's mass itself) that affects the effective local space-time 
structure, see fig.l of ref.[llj. To derive the required correction, let us again denote by {dx, 
dy, dz, dt) the local space-time units of the freely-falling observer S' in the limit M = and 
by 5U the extra Newtonian potential produced by the heavy mass M at the experimental 
set up where one wants to describe light propagation. In a flat-space interpretation, light 
propagation for the S' observer can then be described by the condition 

c^dP 

{ds^)sv = -j^ - {dx^ + df + dz^) = (6) 

where, to first order in 5U , the space-time units (drr, dy, dz, dt) are related to the corre- 
sponding ones (dx, dy, dz, dt) for 6U = through an overall re-scaling factor 

A.:.^ (T, 



and we have also introduced a vacuum refractive index 



A/-=1 + 2M (8) 

Therefore, to this order, light is formally described as in General Relativity where one finds 
the weak-field, isotropic form of the metric 

{ds^)GK = c^dT\l - 2!^) - {dX^ + dY^ + dZ^){l + 2!^) EE c^dr^ - df (9) 

c c 

In Eq.© Un denotes the Newtonian potential and {dT, dX, dY, dZ) arbitrary coordinates 

defined for C/n = 0. Finally, dr and dl denote the elements of proper time and proper length 

in terms of which, in General Relativity, one would again deduce from ds'^ = the same 

universal value c = Jp. This is the basic difference with Eqs.Q-dlj) where the physical 

unit of length is \/ dx'^ + dy'^ + dz"^, the physical unit of time is dt and instead a non-trivial 

refractive index J\f is introduced. For an observer placed on the Earth's surface, its value is 

AA- 1-^^-1.4- 10-9 (10) 



general isotropic metric (A, —B, ~B, —B) depends on two functions which, in a flat-space picture, can 
be interpreted in terms of an overaU re-scaling of the space-time units and of a refractive index. Since physical 
units of time scale as inverse frequencies, and the measured frequencies uj for 5U 7^ are red-shifted when 
compared to the corresponding value uo for 5U — 0, this fixes the value of A. Furthermore, independently of the 
specific underlying mechanisms, the two functions A and B can be related through the general requirement 
AB — 1 which expresses the basic property of light of being, at the same time, a corpuscular and undulatory 
phenomenon [21] . This fixes the value of J\f. 
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where Gn is Newton's constant and M and R are respectively the Earth's mass and radius. 

iii) Differently from General Relativity, in a flat-space interpretation with re-scaled 
units {dx, dy, dz, dt) and M ^ 1, the speed of light in the vacuum c-y no longer coincides 
with the parameter c entering Lorentz transformations. Therefore, as a general consequence 
of Lorentz transformations, an isotropic propagation as in Eq.Q can only be valid for a 
special state of motion of the Earth's laboratory. This provides the operative definition of a 
preferred reference frame S while for a non-zero relative velocity V one expects off diagonal 
elements goi 7^ in the effective metric and a tiny light anisotropy. As shown in Ref.[5], to 
first order in both {J\f — 1) and V/c one finds 



90i 



m - 1 



These off diagonal elements can be imagined as being due to a directional polarization of 
the vacuum induced by the now moving Earth's gravitational field and express the general 
property [22j that any metric, locally, can always be brought into diagonal form by suitable 
rotations and boosts. In this way, by introducing 13 = V/c, k = [M"^ — 1) and the angle 6 
between V and the direction of light propagation, one finds, toO{n) and ©(/J^), the one-way 
velocity [5] 



and a two-way velocity of light 

c,{0) = 



This allows to define the RMS 

Ace ^ 

c 

with 



1 — k/3 cos ( 



J/?2(l + cos2, 



(12) 



c^{9) + Cyin + 9) 



c 

Jf V 



1 



anisotropy parameter B through the relation 

Cy{TT/2 + 9)-Cy{9) 



(C7> 



^^cos(26') 



\B\ 



2 



1 



(13) 



(14) 



(15) 



From the previous analysis, by replacing the value of the refractive index Eq. ()10p and adopt- 
ing, as a rough order of magnitude, the typical value of most cosmic motions V ~ 300 km/s, 
one expects a tiny fractional anisotropy 

Yl 

2 



(Ace) 



\B\ 



0(10-^5) 



c c 

that could finally be detected in the present, precise ether-drift experiments 



(16) 
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3. The present ether-drift experiments 



Let us now consider in more detail the experimental aspects. To increase the statistics, the 
present experiments exhibit rotating optical resonators. For instance, the relative frequency 
shift for the experiment of ref.|25j can be expressed as 

— = S{t) sin 2a;rott + C{t) cos 2u,ott (17) 

where Wrot is the rotation frequency of one resonator with respect to the other which is 
kept fixed in the laboratory and oriented north-south. Notice that, by considering a fully 
symmetric apparatus one should introduce an overall factor of two on the right hand side 
of the above equantion, at least if one wants to express the measured shift in terms of the 
functions S{t) and C{t) extracted from the non-symmetric apparatus of ref. |25] . Notice also 
that in some articles the function S{t) is denoted as B{t). 

In this framework, the existence of possible time modulations of the signal that might 
be synchronous with the Earth's rotation has always represented a crucial ingredient for 
the analysis of the data. This expectation derives from a model where one assumes a fixed 
preferred frame S. Then, for short-time observations of 1-2 days, the time dependence of a 
hypothetical physical signal can only be due to (the variations of the projection of the Earth's 
velocity V in the interferometer's plane caused by) the Earth's rotation. In this case, the two 
functions S{t) and C(t) admit the simplest Fourier expansion [25] {t' = Usidt is the sidereal 
time of the observation in degrees) 

S{t) = So + Ssi sin t' + Sci cos t' + Ss2 sin(2t') + Sc2 cos(2t') (18) 

C{t) = Co + Csi sin t' + Cci cos t' + Cs2 sin(2t') + Cc2 cos(2t') (19) 

with time-independent and Fourier coefficients. 

This theoretical framework, accepted so far by all experimental groups, leads to average 
the various Ck and Sk obtained from fits performed during a 1-2 day observation period. 
By further averaging over many short-period experimental sessions, the data support the 
general conclusion [Ml [271 [28] that, although the typical instantaneous S{t) and C{t) are 
indeed 0(10"^^), the global averages (Cfc)^^^ and {SjS)^^^ for the Fourier coefficients are 
much smaller, at the level 0(10^^^), and, with them, the derived parameters entering the 
phenomenological SME |29| and RMS models. 

However, there might be different types of ether-drift where the straightforward parame- 
terizations Eqs. (fT8]l . ([Hj) and the associated averaging procedures are not allowed. In fact, 
before assuming any definite theoretical scenario, one should first ask: if light were really 
propagating in a physical medium, an ether, and not in a trivial empty vacuum, how should 
the motion of (or in) this medium be described? Namely, could this relative motion exhibit 
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variations that are not only due to known effects as the Earth's rotation and orbital rev- 
olution? As discussed in [5], by representing the physical vacuum as a fluid, the standard 
assumption of smooth sinusoidal variations of the signal, associated with the Earth's rotation 
(and its orbital revolution), corresponds to assume the conditions of a pure laminar flow 
associated with simple regular motions. Instead, by adopting the model of an underlying 
turbulent medium there might be other forms of time modulations. In this alternative sce- 
nario, the same basic experimental data might admit a different interpretation and a deflnite 
instantaneous signal Ai^{t) ^ could become consistent with (Ck)^^^ ~ (SkY^^ ~ 0. 
To discuss this alternative scenario, it is convenient to first re- write Eq. (|17p as 

^ = A{t) cos{2ujrott - 2eo{t)) (20) 

where 

C{t) = A{t) cos 2eo{t) S{t) = A{t) sin 29o{t) (21) 

so that 



Ait) = ^S^t) + C^t) (22) 

Here 9o{t) represents the instantaneous direction of a hypothetical ether-drift effect in the x-y 
plane of the interferometer (counted by convention from North through East so that North 
is ^0 = and East is = ^/2)- By also introducing the magnitude v = v{t) of the projection 
of the full V, such that 

Vx{t) = v{t) sin 0o{t) Vyit) = vit) cos 6'o(t) (23) 

we obtain the theoretical relations [5] 

A{t) = l\B\^-^ (24) 
2 



and 

2 c 2 c 



In the forthcoming section we shall produce a numerical simulation by assuming a model of 
turbulent flow for the velocity components Vx{t) and Vy{t). 



4. Numerical simulation of a physical, stochastic component 

Before trying to simulate a physical stochastic component of the signal, to obtain the correct 
normalization, we should first subtract from the existing data the known spurious effects. To 
obtain a precise statistical indicator we shall consider the root square of the Allan variance 
(RAV) for an integration time r ~ 1 second which we'll take as our definition of instantaneous 
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signal. In fact, our model predicts typical frequency shifts of about 1 Hz so that it only makes 
sense to compare with sequences of data collected at time steps of 1 second or larger. The 
RAV describes the time dependence of an arbitrary function z = z{t) which can be sampled 
over time intervals of length r. In this case, by defining 

z{ti-T) = - r dt z{t) = Zi (26) 
^ Ju 

one generates a r— dependent distribution of values. In a large time interval A = Mr, the 
RAV is then defined as 

RAV(t) = ^/a^{z,T) (27) 

where 

M 

a\z,T) = —Y,{z,--^^+lf (28) 
1=1 

Now, for the non-rotating set up, the RAV of the frequency shift for r ~ 1 second was 
determined [9] to be 0.8 Hz (2.8 • 10~^^ in dimensionless units) and found much larger 
than the corresponding disturbances in the individual resonators (typically about 0.02-0.03 
Hz). The only exception is the possible effect of thermal noise in the mirrors of the optical 
resonators. This particular component, which should be independent of the integration time, 
on the basis of the results of ref.pO], was estimated in ref.|31j to be about 1.15 • 10^^^ in 
dimensionless units. Therefore, for a laser frequency = 2.82 • 10^^ Hz, we would expect 
RAV(thermal-noise)~ 0.32 Hz. It is questionable how to subtract this effect from the full 
measured value 0.8 Hz. One might argue that, if the physical signal has also a stochastic 
nature, one should subtract in quadrature. This would give 

RAV(physical, r ~ 1 second) = ^(0.8)2 - (0.32)2 _ o.73 Hz (29) 

Instead, we shall adopt the more conservative attitude of subtracting linearly, i.e. 

RAV(physical, r ~ 1 second) = 0.8 Hz - 0.32 Hz = 0.48 Hz (30) 

or 1.7 • 10~^^ in dimensionless units. Since for a symmetric non-rotating set-up the physical 
signal is simply 2C(t)t'o, we conclude that there is a potentially important contribution to 
C{t) which corresponds to a stochastic signal with an Allan variance of about 8.5 • 10~^^ for 
r ~ 1 second. This value will be the basic input for our simulation. 

Let us now return to Eqs. (p5|) and assume for the velocity components Vx{t) and Vy{t) 
a toy model of turbulent flow. To this end, we shall follow ref.jSSj where velocity flows, in 

''We tried to obtain an analogous indication from the other experiment of ref |27). However, for r ~ 1 
second, it is not so easy to determine the value of the RAV. In fact, by inspection of their figure 2, in the 
narrow range from r = 0.8 seconds to r = 1 second, the data for the non-rotating set-up (the red dots) exhibit 
a very steep, sizeable decrease from about 2.8 ■ 10~'^ down to 1.4 ■ lO"'^^. 
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statistically isotropic and homogeneous 3-dimensional turbulence, are generated by unsteady 
random Fourier series. The perspective is that of an observer which moves in the turbulent 
fluid and wants to simulate the two components of the velocity in his x-y plane at a given 
fixed location in his laboratory. This leads to the general expressions 

oo 

= '^[Xn{l) COSOJnt + Xn{2) s'mUJnt] (31) 
n=l 

oo 

^2/(*) = [yn(l)cosa;„t + y„(2)sina;„t] (32) 

n=l 

where ujn = 2n7r/T, T being a time scale which represents a common period of all stochastic 
components. In our simulation we have fixed the typical value T = ?day= 24 hours. However, 
we have also checked with a few runs that the statistical distributions of the various quantities 
do not change appreciably if we vary T in the rather wide range 0.1 Tjay < T < 10 T^ay 

The coefficients Xn{i = 1,2) and y„(i = 1,2) are random variables with zero mean. They 
have the physical dimension of a velocity and we shall denote by [—v, v] the relevant interval 
of these parameters. In terms of v the quadratic mean values can be expressed as 

{xl{i = 1, 2)) = {yl{i = 1, 2)) = ^ (33) 

for the uniform probability model (within the interval [— v,^]) which we have chosen for 
our simulations. Finally, the exponent /3 controls the power spectrum of the fluctuating 
components. For our simulation, between the two values /3 = 5/6 and /? = 1 reported in 
ref.[32], we have chosen (3 = 1 which corresponds to the point of view of an observer moving 
in the fluid. 

Thus, within this simple model for the stochastic signal, v is our only free parameter 
and will be fixed by imposing that the generated C— values give a RAV of 8.5 • 10~^^ for 
integration time r = 1 second. By taking into account the typical variation of the results, 
due to both the truncation of the Fourier modes and the dependence on the random sequence, 
this constraint gives a range v ~ (332 ± 8) km/s which, remarkably, has a definite counter 
part in the known Earth's cosmic motion with respect to the Cosmic Microwave Background 
(CMB). In fact, it coincides exactly with the daily average of the projection (u^) ~ 332 
km/s in the interferometer's plane for an apparatus at the latitude of Berlin-Diisseldorf. This 
can be checked by using the relation [5] 

{v"^) = ^1 — sin^ 7 cos^ ^ ~ 2 ^ ^ ^"^^^ 

and setting V = 370 km/s, angular declination 7 ~ —6 degrees and co-latitude x ~ 38 
degrees. 



9 



50130 



50140 



50150 



50160 



S(t) 0.. 



I I \ I \ I \ 

50130 50140 50150 50160 

time (s) 



Figure 1: A simulation of the instantaneous values of the C and S functions, in units 10 , 
as obtained from a typical sequence of 40 seconds . 

After these preliminaries, the results of our numerical simulation can be illustrated by 
starting from the building blocks of our scheme, namely the instantaneous values Cj = C(ti) 
and Si = S(ti) of the C and 5— functions. In Fig.l we report a typical sequence of 40 values 
of these functions. 

In terms of these basic quantities, one can construct a first type of averages over a time 
scale T = N seconds 

i+N-l i+N-l 

C{tv,N) = - Cn S{tf,N) = - Sn (35) 

n=i n—i 

so that C{ti; 1) = Cj and S{ti\ 1) = Si. This first type of averaging is essential to compare with 
the experimental data where the C and S*— functions are always determined after averaging 
the basic instantaneous data over times r in the typical range 40-400 seconds. With these 
auxiliary quantities, collected during a large time scale A = Mr, one can form a statistical 
distribution and determine mean values 

M M 

{C)r = j^Y.'^{ti-T) ^S)r = jjY^iti-r) (36) 



i=i 1=1 



and variances 



^ {C{tr,T)-{C)r)' ^ {Sit,;r)-{S)ry 

-c{r) = Y M^I -sir) = Y 

i=i 1=1 



(37) 
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Figure 2: We show, see (a) and (h), the histograms of the simulated C and S values, in 
units lO""*^^, for r = 1 second. The vertical normalization is to a unit area. The mean values 
are {C)r = -1.1 • 10"^^, {S)r = -1.9 • 10^^*^ and the standard deviations adr) = 8.5 • lO"^*^, 
^^si^) = 9.4 • 10^^^. We also show, see (c) and (d), the corresponding plots in logarithmic 
scale and the fits with Eg. Ii3^). The total statistics correspond to a time A = Mr =864-00 
seconds. 
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Figure 3: The histograms of the simulated C and S values, in units 10""^^, and the cor- 
responding gaussian fits for r = 40 seconds. The vertical normalization is to a unit area. 
The mean values are {C)r = —9 • 10~^^, {S)r = —5 • 10~^^ and the standard deviations 
(^c{t) = 1-34 • 10~^^, ^"^(t) = 1.48 • 10""*^^. The total statistics correspond to a time 
A = Mr =864000 seconds. 

We report in Fig. 2, for r = 1 second, the distribution functions of the simulated C and 
S values (panels (a) and (b)). Notice that these distributions are clearly "fat-tailed" and 
very different from a Gaussian shape. This kind of behavior is characteristic of probability 
distributions in turbulent flow at small time scales (see e.g. |33tl34)). To better appreciate the 
deviation from Gaussian behavior, in panels (c) and (d) we plot the same data in a log— log 
scale. The resulting distributions are well fitted by the so-called (7— exponential function j35j 

y = ^(1 - (1 - ^)x5)i/(i-5) (38) 

with an entropic index q= 1.01. 

By starting to average the instantaneous values, the statistical distributions of the simu- 
lated C and S tend to assume a gaussian shape. This is already evident from about r = 5 — 6 
seconds. In Fig. 3 we show the two distributions for r = 40 seconds. 

As it might be expected, for all r the statistical averages {C)r and {S)r are vanishingly 
small in units of the typical instantaneous signal 0(10^^^) and any non-zero value has to 
be considered as statistical fluctuation. The standard deviations, on the other hand, have 
definite values and exhibit a clear 1/\/t trend so that, to good approximation, one can express 

, , 8.5-10-16 , , 9.4-10-16 , , 

^c(r) ■==- as{T) ■==- (39) 

Vr(sec) V^(sec) 

By keeping v fixed at 332 km/s, the above two values for r = 1 second have an uncertainty 
of about 5% which reflects the typical variation of the results due to both the truncation of 
the Fourier modes and the dependence on the random sequence. 



12 



Notice that our model predicts a monotonic decrease of the dispersion of the data by 
increasing the averaging time r and, therefore, cannot reproduce the linear increase of the 
Allan variance which is seen, in all present room temperature experiments, above about 
T = 100 seconds. However, this might also be a spurious effect. For instance, in the cryogenic 
experiment of ref . |36j the Allan variance (in the quiet phase between two refillings of the tank 
of liquid helium) was found to exhibit a monotonic decrease up to about r = 250 seconds. It 
remains to be seen if this decreasing trend will be confirmed by the forthcoming generation 
of experiments with cryogenic sapphire resonators |17j that are expected to have a short-time 
stability of a few 10~^^. Thus it will be possible to obtain a precise check of our predictions. 
In particular, the typical instantaneous signal should be about 200—300 times larger than 
the experimental sensitivity and the distributions of the C(ti;r) and S'(ti;r), for r = 100 
seconds, should extend up to values which are still 20—30 times larger. 

5. Summary and conclusions 

As discussed in the Introduction, ether-drift experiments play a fundamental role for our 
understanding of relativity. In fact, they represent the only known experiments which, in 
principle, can distinguish Einstein's interpretation from the Lorentzian point of view with a 
preferred reference frame S. Up to now, the interpretation of the data has been based on a 
theoretical model where the hypothetical S is assumed to define a fixed direction in space. 
In this way all type of signals that are not synchronous with the Earth's rotation tend to be 
considered as spurious instrumental noise and no particular effort is made to understand if 
there could be genuine physical effects which do not fit within the adopted scheme. 

However, there is a logical gap which has been missed so far. Even though the relevant 
Earth's cosmic motion corresponds to that indicated by the anisotropy of the CMB (V ~370 
km/s, angular declination 7 ~ —6 degrees, and right ascension a ~ 168 degrees) it might 
be difficult to detect these parameters in microscopic measurements of the speed of light 
performed in a laboratory. The link between the two concepts depends on the adopted model 
for the vacuum. The point of view adopted so far corresponds to consider the vacuum as 
some kind of fluid in a state of regular, laminar flow. In these conditions global and local 
informations on the flow coincide. 

We believe that without fully understanding the nature of that substratum that we call 
physical vacuum, one should instead keep a more open mind. For instance the physical 
vacuum might be similar to a form of turbulent ether, an idea which is deep rooted in the 
basic foundational aspects of both quantum theory and relativity (see e.g. [llj) and finds 
additional motivations in those representations of the vacuum as a form of 'space-time foam' 
which indeed resembles a turbulent fluid [12l [131 El IISl US]- In this case, global and local 
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indications might be very different and there could be forms of random signals that have 
a genuine physical origin. For instance, by combining the point of view of ref.[5j, where 
gravity is considered a long-wavelength phenomenon which emerges from a space-time which 
is fundamentally flat at very short distances, with the idea of a turbulent ether, one arrives 
to a stochastic signal of typical magnitude 10~^^ which could fit very well with the present 
experimental data. 

For this reason, after reviewing in Sects. 2 and 3, the general theoretical framework and 
the basics of the present ether-drift experiments, we have presented in Sect. 4 a simple nu- 
merical simulation of the possible effects that one might expect in a statistically isotropic 
and homogeneous turbulent flow. After subtracting the known forms of disturbances, we 
have found that the observed distribution of the data over a time scale of 1 second require 
a value ~ 332 km/s of the scalar velocity parameter which characterizes the fluctuations. 
Remarkably, this has a definite counter part in the known Earth's cosmic motion with re- 
spect to the CMB. In fact, it corresponds exactly to the daily average of the projection of 
the Earth's velocity in the interferometer's plane for an apparatus placed at the latitude of 
Berlin-Diisseldorf. However, this correspondence with the global Earth's motion is only valid 
at the level of statistical distributions and is not detectable from the naive time dependence 
of the data. In conclusion, we predict a typical instantaneous signal which should be about 
200—300 times larger that the sensitivity of the forthcoming generation of cryogenic exper- 
iments |17) . A confirmation of this prediction would represent unambiguous evidence for a 
physical ether-drift with non-trivial implications for our understanding of both gravity and 
relativity. 
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